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We report planar tunneling measurements on thin films of YBa2Cu307_a; at various doping levels 
under magnetic fields. By choosing a special setup configuration, we have probed a field induced 
energy scale that dominates in the vicinity of a node of the d-wave superconducting order parameter. 
We found a high doping sensitivity for this energy scale. At Optimum doping this energy scale is 
in agreement with an induced id xy order parameter. We found that it can be followed down to low 
fields at optimum doping, but not away from it. 



PACS numbers: 74.72.Bk, 74.50.+r 
I. INTRODUCTION 

It is by now well established that cuprate supercon- 
ductors have an order-parameter with a dominant d- 
wave symmetry, characterized by node-lines along the 
[110] and equivalent directions 1 -. Nodal quasi-particles 
become then the dominant low energy excitations^. In- 
teresting phenomena have been predicted to occur when 
the d-wave superconductor is subjected to a magnetic 
field perpendicular to the superconducting planes: an 
energy gap should develop at the nodes, that increases 
with the square root of the applied field. This has been 
explained by Laughlin 3 who assumes an additional imag- 
inary id X y component which increases with the magnetic 
field. Discrete Landau energy levels were also predicted 
to develop in the nodal regions by Gor'kov-Schrieffer 4 
and Anderson^. However, the observation of nodal finite 
energy levels has encountered theoretical and experimen- 
tal difficulties. 

It has been argued, that in the mixed state, superfiuid 
screening currents result in a Doppler-shift of the Lan- 
dau levels larger than the level spacing, rendering their 
observation impossible^^. This Doppler shift will also 
obscure a possible id xy component^. Tunneling experi- 
ments performed along a nodal direction, which should 
in principle be an ideal method to probe nodal states, are 
dominated by low energy Andreev - Saint James (ASJ) 
surface states due to the d-w&ve symmetry, resulting in a 
characteristic Zero Bias Conductance Peak (ZBCP). The 
degeneracy of the Andreev - Saint-James states is lifted 
by screening currents splitting the ZBCP, an effect that 
can be confused with that of nodal finite energy levels. 
Spontaneous time reversal symmetry breaking effects are 
also sometimes observe d 10 ' 11 : 12 ! 13 ! 14 . In addition, it is not 
trivial to distinguish between the predictions of the Lan- 
dau States and of the minority order parameter theories. 
This is because, the energy of the first Landau level is 
equal to the amplitude of the id xy component predicted 
by Laughlin. 

To address these difficulties we have used field cycles 
that enable us to distinguish between the Doppler shift 



and other possible spectral contributions. We have con- 
cluded that data taken in decreasing fields is essentially 
free of Doppler shift effects and can be used to identify 
finite energy nodal states. We confirmed that this energy 
follows the predicted square root of field behavior at op- 
timum doping. By extending our measurements to 20 
mK, we were able to follow the evolution of these states 
down to fields of the order of a few 1000 Gauss, where 
the Landau level interpretation is excluded due to the 
long length of the corresponding trajectories. Finally we 
report that the doping level had a strong influence on the 
splitting behavior of the ZBCP. The square root behavior 
is not obeyed at low fields when deviating from optimum 
doping. Such deviations were reported befor o 14 : 15 , how- 
ever, the data was taken in increasing fields and was an 
admixture of the studied phenomena and Doppler shift. 
In the new data presented and studied in this paper, the 
Doppler shift contribution is essentially eliminated, and 
the doping effect appears far more clearly. In Summary, 
our results favor the existence of additional id xy order 
parameter component predicted by Laughlin rather than 
that of the formation of nodal energy levels predicted by 
Gor'kov-Schrieffer— and Anderson^. 

This paper will be organized as follows: we begin with 
a theoretical background of the Doppler-shift effect. In 
Sec. IHII we present our experimental setup enabling us to 
distinguish between the two contributions. Our tunneling 
results at optimum doping will be shown together with 
the low temperatures measurements at 20mK (Sec. II VI) . 
We then compare our results to theory in Sec. [V] and 
finish with our conclusions and findings. 

II. THEORETICAL BACKGROUND 

We wish to discuss the differences between two the- 
oretical approaches regarding the development of finite 
nodal energy states under applied magnetic fields. 

In the first approach, by Laughlin 3 -, the free energy 
of a rf-wave superconductor subjected to a magnetic field 
perpendicular to the superconducting planes can be min- 
imized by the inclusion of an additional id xy component 
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FIG. 1: (color online) Schematic illustration of the electronic 
momentum space in a d-wave superconductor under an ap- 
plied magnetic field. The x and y axis are parallel to the 
[100] and [010] crystallographic directions perspectively. For 
simplicity we assume: a cylindrical Fermi surface, a d x 2_ y 2- 
wave superconducting gap and nodes at ±45° from the prin- 
cipal axes, and a magnetic field, B, parallel to the z axis. The 
quasi-particles cycle of multiple Andreev - Saint- James reflec- 
tions forming the nodal energy level is marked by the solid 
line. As an alternative theoretical description for the energy 
scale an induced id xy order parameter (marked in purple) 
is predicted to develop mainly in the vicinity of the nodes. 
The collimation of the injected electrons in a planar tunnel- 
ing configuration of the experiment is marked by the green 
triangles for two different junction orientations. The field in- 
duced nodal energy scale can only be probed for tunneling 
along the node direction (left triangle). 



to the main d x 2_ y ? component (illustrated in Fig. [T]by 
the id xy order parameter marked in purple). The id xy 
component breaks the symmetry in such a way that op- 
posite currents will flow on opposite faces of the sample 
creating a magnetic moment parallel to the applied field. 
If the moment and the applied field are parallel to each 
other the free energy will be minimized. 

In the second approach, following Anderson^, we con- 
sider the motion of a quasi-particle in a nodal region (Fig. 
[I}. Under an applied magnetic field B, it acquires a ve- 
locity component parallel to the Fermi surface. In the 
absence of a superconducting order parameter, it will be 
in one of the Landau levels, determined by the cyclotron 
frequency uj c = -Nr-, where e is the quasi-particle charge, 
c is the speed of light and m* is the effective electron 
mass. However, when hio c < A, where A is the ampli- 
tude of the superconducting gap, the usual cyclotron mo- 
tion is not possible. Instead, as schematically described 
in Fig. [T] a series of Andreev - Saint- James reflections 
in momentum space will occur—. This process can only 
occur at certain energy levels for which the total phase 
change during a Saint-James cycle is a multiple of 2tt. 
These energy levels correspond values of the current that 
flows around the Fermi surface. 

Interestingly, the amplitude of the minority component 
is the same as the energy of the first Landau level 4 ^. The 



two approaches lead to the same energy exactly: 

e(B) = ±2y/hw c A (1) 

Nonetheless, there are substantial differences between 
these two approaches. While the Gorko'v-Schrieffer- 
Anderson theory assumes that the order parameter is not 
altered by the magnetic field, its modification is a key 
prediction in Laughlin's theory. The latter also predicts 
a transition temperature above which the <i-wave sym- 
metry is recovered. Finally Gor'kov-Schrieffer-Anderson 
predict a series of energy levels while Laughlin only pre- 
dicts a finite gap value. 

Tunneling along a nodal direction of a d-wave super- 
conductor is done through a surface where zero energy 
bound states are present due to the interference of quasi- 
particles that undergo Andreev - Saint-James reflections 
from lobes of the order parameter having phases that 
differ by 7 r 13 i 16 . These bound-states should appear as a 
conductance peak at zero bias in an in-plane tunneling 
spectrum^. As shown by Fogelstrom et ai, this zero bias 
peak splits into two spectral peaks due to a Doppler-shift 
from superfluid currents flowing parallel to the surface 1 ^. 
The peaks bias are proportional to v s ■ p F , where v s is 
the superfluid velocity and p F the Fermi momentum of 
the probed states. For example, when a magnetic field, 
H, is applied parallel to the surface, Meissner currents 
Doppler-shift will produce spectral peaks which are lin- 
ear with H up to a field of the order of the thermody- 
namical critical field where saturation is reached (about 
1 Tesla in the case of YBaaCuaOr-z (YBCO)). 

Since a nodal energy scale and a Doppler shift will 
both split the zero bias conductance peak under an ap- 
plied field, as has been observe d 10 ' 11 ^ 2 ' 14 , one must find a 
method to distinguish between both mechanisms and de- 
termine the difference in the predicted field dependences. 
An obvious difference between them is that field induced 
nodal energy scales are best observed in the absence of su- 
perfluid currents, while a Doppler shift effect exists only 
in their presence. 

A method, which we have already usecU^, consists in 
performing magnetic field cycles. Meissner currents are 
quite different in increasing and decreasing fields because 
the Bean-Livingston barrier which can retard the pene- 
tration of vortices through strong surface currents, up 
to a field of the order of the thermodynamical critical 
field, is only effective against flux penetration (increas- 
ing fields) and not against flux exit (decreasing fields). 
As a result, strong surface Meissner currents, on the 
scale of the London penetration depth, exist only in in- 
creasing field o 19 ' 20 i 21 . Other types of current that can 
produce a Doppler shift are screening currents around 
vortices^ 2 - and Bean's critical state currents 23 . The lat- 
ter reversesign with field reversal, and in high fields, ex- 
tends into the entire thickness of the sample. They are 
typically weaker than the Meissner currents in the Bean- 
Livingston regime. 

In planar tunneling experiments, electrons are injected 
across a dielectric barrier into the superconductor. The 
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transmission probability decays exponentially with the 
increasing angle between the electron's momentum and 
the normal to the interface resulting in a collimated cur- 
rent. In a typical junction the momentum divergence has 
a angle of 10-20 degrees known as the tunneling cone. 
The width of this cone, which can vary slightly from one 
junction to another, will influence the Doppler shift of 
the zero energy surface bound states, however, it will 
not modify the energy of nodal states nor that of an in- 
duced id xy order parameter component. 

In summary, a zero bias conductance peak is expected 
to appear in a d-wavc tunneling along the node direction 
and to split into two spectral peaks when a magnetic field 
is applied perpendicular to the CuC>2 planes due to a field 
induced nodal energy scale or via a Doppler-shift effect. 
The two effects behave differently in a magnetic field. In 
the next section, we show a way to minimize the Doppler- 
effect which allowed us to probe the field induced nodal 
energy scale alone. 



III. EXPERIMENTAL 

Thin YBCO films were grown using DC off-axis sput- 
tering deposition. In order to minimize (103) oriented 
grains a buffer layer of PrBa2Cu307_ :r was first deposited 
using RF off-axis sputtering on top of the substrate 2 ^. 
We used a SrTiC-3 and LaSrGaCU substrates for the (110) 
and (100) oriented films respectively. These films have a 
well defined [001] direction parallel to the surface of the 
film. 6 — 29 x-ray diffraction patterns showed the rele- 
vant peaks for the desired orientation. Scanning electron 
microscopy and atomic force microscopy showed a well 
defined crystallographic growth and surface roughness of 
a few nanometers. Resistivity measurements showed the 
expected in-plane anisotropy. In addition, the tempera- 
ture dependence of the ab plane resistivity allowed us to 
estimate the doping level in our films. It changes from a 
positive curvature for overdoped, linear with temperature 
for optimally doped and negative curvature for under- 
doped film o 25 ' 26 ! 27 ' 28 . I-V characteristics were measured 
using a current source and a digital voltmeter. The tun- 
neling conductance spectra were calculated by differenti- 
ating the I(V) curves. Table U shows the characterization 
values for the 15 representative junctions used in the fig- 
ures of this paper. 

In our experiments the tunneling junction is created 
by placing an indium electrode on top of the surface 
of freshly prepared thin YBCO films. At the metal- 
superconductor interface a thin insulating indium-oxide 
layer is then formed, which is stable over weeks and many 
thermal cycles. We can verify the quality of the tunneling 
junction by several methods. First, by lowering the tem- 
perature below the indium critical temperature we can 
measure the well-know indium tunneling spectrum which 
dominates the low energy spectra. Also the indium spec- 
trum disappears as we increase the magnetic field or heat 
up the sample above the indium's critical field and tem- 



perature, respectively. Second, we ensure that the high 
bias conductance is insensitive to the magnetic field and 
temperature below the YBCO critical temperatureas, as 
expected for tunneling spectroscopy. 

Our sample and the field configuration are favorable 
for several reasons. First, the magnetic field is applied 
parallel to the surface (which avoids threading the tun- 
nel junction with vortices) and at the same time per- 
pendicular to the Cu02 planes, as required for the ob- 
servation of finite energy nodal states. Second, the un- 
desired Doppler-shift due to superfluid currents is mini- 
mized because in our geometry v s ■ p F is at a minimum 
since the dominant currents flow parallel to the inter- 
face. Third, by comparing data taken in increasing and 
decreasing fields it is possible to identify the effect of 
the Doppler-shift due to strong surface Meissner cur- 
rents that exist only in increasing field a 19 i 20 ' 21 . There- 
fore, measurements in decreasing fields are less affected 
by the Doppler-effect. Also, the intensity of the Meiss- 
ner currents themselves can be minimized by using films 
whose thickness is smaller than the London penetration 
dept h 18 ' 29 . Finally, we have performed experiments both 
on (110) and (100) orientated films. While for the [110] 
direction we expect to probe the field induced energy 
scale, for the [100] direction it should not be observed. 
We emphasize that in our method due to the planar ge- 
ometry of the junction and the resulting tunneling cone, 
we probe a specific direction in fc-space rather than the 
k averaged local density of states as probed by scanning 
tunneling microscopy. 



IV. RESULTS 
A. Probing the field induced nodal energy scale 

In Fig. [5] we show tunneling spectra for an optimally 
doped (110) oriented film as a function of magnetic field 
applied perpendicular to the Cu02 planes. We notice 
that the spectral peak bias values, S(H), are always larger 
in increasing rather than in decreasing fields. In addition, 
the peak seen in decreasing fields is well defined for all 
fields, while that seen in increasing fields becomes very 
broad and is too broad to be identified (Fig. \2§ , while in 
decreasing fields it remains well defined up to more than 
22 T (see for example Fig. [7]). 

We demonstrate in Fig. [3j and [4] that the spectral 
peak bias value does not increase linearly at low fields, 
and does not saturate at high field. This is in contradic- 
tion with the Doppler shift theory of Zero Energy Sur- 
face Bound States. A substantial difference between the 
behavior of (100) and (110) oriented films is observed. 
While for thin (110) films, the spectral peaks are clearly 
seen even at low fields, no such peaks can be detected in 
thin (100) films. Since in (100) films, the only splitting 
mechanism is the Doppler shift effect, the absence of the 
spectral peaks in such films is indicative of the insignif- 
icance of that effect in films thinner than the London 
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TABLE I: Samples characterization. As explained in the text, from the resistivity temperature dependence measurement, R(T), 
we can estimated the samples doping. The zero field spectral peak bias value is noted as So. T c is the temperature at zero 
resistivity. The transition temperature width is determined by a Gaussian fit to ^ at the transition. 
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penetration depth. Upon increasing the (100) oriented 
films thickness, the contribution of the Doppler shift is 
also increased and the spectral peaks are recovered^. 

The third qualitative evidence supporting the ar- 
gument that we probe a field induced nodal energy 
scale rather than a Doppler-shifted zero bias conduc- 
tance peak comes from the shape of the tunneling con- 
ductance at zero bias. While the Doppler shifted ZBCP 
results in a V-shape conductance^^, one expects a U- 
shaped conductance at zero bias for the field induced 
nodal energy scale — . Because it is difficult to distin- 
guish between these two shapes due to thermal popula- 
tion effects, one should perform measurements at very 
low temperatures. In Fig[4] we show tunneling spectra 
taken at 20mK under high magnetic fields. The ob- 
served U-shaped conductance is in agreement with the 
field induced nodal energy scale scenario and contrasts 
the Doppler shift model. 

We have therefore demonstrated that a Doppler shift 
due to Meissner screening currents does not play an im- 
portant role in our tunneling measurements of thin films 
in decreasing fields. However, the effect of a nearby 
vortex^, or Bean's critical state currents^ could still 
take place in a decreasing magnetic field. These currents 
reverse their polarity when decreasing the magnetic field 
from the maximum value reached. If they dominate, the 
total current should be zero at some field. In this case, a 
zero bias peak should reappear at a finite magnetic field. 
Such a behavior was never observed, ruling out that the 
spectral peaks in decreasing fields arise from the effect of 
a nearby vortex or a strong Bean critical state currents 
(see for example a field cycle in Fig. [5]). 

Finally we note that the peak position measured in 
decreasing fields is extremely reproducible, for a variety 
of samples and junctions. This is in contrast with the 
expected behavior of the Bean critical currents and the 
Doppler shift effect that are strongly dependent on film 



thickness, tunneling cone and surface barrier formation. 

In contrast with the difficulties encountered in trying 
to explain the data shown in Fig. [2] and Fig. [3] by a 
Doppler shift of zero-energy surface bound-states, a field 
induced nodal energy scale provides a reasonable expla- 
nation. As shown in Fig. [H the data taken in decreas- 
ing fields fits the predicted square root dependence on 
the magnetic field. The field hysteresis shown in Fig. [5] 
and Fig. [3] can be understood as due to a Doppler shift 
by Meissner currents in increasing fields; the peaks are 
shifted to higher bias and are broadened until they can- 
not be identified anymore in very high fields. 

To summarize, we have shown that for an optimally- 
doped sample, tunneling measurements on (110) oriented 
film reveal an energy scale that can be understood cither 
as an additional complex order parameter in the form of 
id xy or as the first Landau state in the vicinity of the 
d-wave node. 



B. Implausibility of nodal Landau levels 

After ruling out the Doppler shift effect, we now con- 
centrate on distinguishing between the nodal Landau 
states approach of Gor'kov-Schrieffer— and Anderson^ 
and the nodal order parameter component predicted by 
Laughlin^. Although the first Landau level and the nodal 
order parameter have identical field dependences, the fol- 
lowing evidences favor the later. 

An important feature is observed in Fig. [4] only two 
peaks at ±<5 are visible. This is incompatible with the 
theories in refs4^, which predict a series of energy levels 
that should manifest themselves as peaks in the tunnel- 
ing spectrum. These peaks have never been observed. 
One could argue that the higher energy levels are hid- 
den due to scattering. The visibility of the first peak 
at relatively high temperatures (T > A.2K) under small 
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FIG. 2: (color online) Tunneling spectra, dI/dV(V), for var- 
ious magnetic fields (sample SI, at 0.3 K). Increasing from 
zero magnetic field (black) and decreasing from 16 T (red), 
they present different tunneling spectra in a given magnetic 
field. We define 28 as the distance between the two maxima 
in the spectrum. For a given magnetic field, 8(H) is always 
larger in increasing fields than in decreasing one. This is due 
to the Doppler-shift effect resulting in shifting and smearing 
the field induced nodal energy scale peaks to higher biases 
in increasing magnetic fields. In fields higher than 4 T, the 
peaks can only be identified in decreasing magnetic fields. 



magnetic fields (H ~ 0.3T) suggests that the scattering 
processes are not strong enough to obscure the higher 
peaks at 20mK and 9T, yet, these peaks are not observed 
in Fig. SI 

Further evidence comes from estimating the trajectory 
length at the node area. If we define 9 as the angle from 
the anti-nodal direction at which the particles meet the 
superconducting gap A and perform an Andreev - Saint- 
James reflection, and ±a as the trajectory angle mea- 
sured from the node we get a — 2(45 — 9). Using Eq. Q] 

we calculate a for small angles to be a = = 2^^^. 
In the node area, the trajectory can be treated classically 
using the cyclotron frequency oj c . In the node vicinity, 
the trajectory time t = — gives the trajectory length x: 



, kmc T 



(2) 



where Vf is the Fermi velocity. Using Eq. [21 we calcu- 
late that for B = 1 Tesla, A = 20meU and Vt = 10 6 ^- 

' J sec 

the trajectory length is 8600A which is several times 
larger than the film's thickness. Scattering from the sur- 
face will not allow the formation of Landau states, unless 
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FIG. 3: (color online) Spectral peak bias value, 8, as a 
function of magnetic field at 4.2 K. Black triangles (dia- 
monds) represent the peak positions in decreasing (increas- 
ing) magnetic fields for a (110) oriented film having thickness 
of 1600A(sample S2). While a Doppler-shift m increasing 
fields prevents observation of the field induced nodal energy 
scale, the measurements in decreasing fields are free of the 
Doppler-shift effect and allow unambiguously identification 
of the field induced nodal energy scale. Dashed purple line is 
a fit to 8 = AH*, with A = 1.02± 0.05 meV/T5 in excel- 
lent agreement with Eq.(l). The red circles represent (100) 
oriented film with a thickness of 1000A(sample S3). The field 
induced nodal energy scale are not observed, as expected. 



the applied field is of the order of 100T. But in fact the 
nodal scale is observed at low temperatures (20mK) down 
to a fraction of a Tesla (Fig0|). 

Two additional indications favoring Laughlin's theory 
over Gor'kov-Schrieffer- Anderson will be discussed in the 
the following sections: the effects of doping in section 
II V CI and a first order phase transition of the field induced 
nodal energy scale measured by Elhalel et al^ in section 



Effect of doping on field induced nodal energy 
scale 



1. Underdoped case 

For underdoped samples at zero field only a single 
peak is observed at zero bias (the ZBCP). We shall now 
demonstrate that this is not due to the thermal smear- 
ing of two peaks at finite energy. To do that we reduced 
the temperature to 0.3K while applying a small mag- 
netic field perpendicular to the c-axis (and parallel to 
the sample's surface). This field quenches superconduc- 
tivity in the indium counter electrode but has no effect 
on the ZBCP as has been demonstrated by Krupke and 
Deutscher— . The results are shown in Fig. [5] The sam- 
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FIG. 4: (color online) (a) Tunneling spectra taken at 20mK 
in decreasing magnetic fields (sample S4). Note the constant 
conductance in the vicinity of zero bias supporting the exis- 
tence of an energy scale e = ±8(H) (b) Spectral peak value in 
decreasing fields versus square root of the applied magnetic 
field measured at 20mK. the line is a fit to high fields having 
a slope of 1.01 meV/T^. 
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FIG. 6: Spectral peak value, 5, as a function of the square- 
root applied magnetic field (sample S6). The line is a linear fit 
to the high magnetic field data. It has a slope of 0.9 meV/T^ . 
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FIG. 5: Tunneling conductance for a (110) film taken at 0.3K 
(sample SI). The black line measured in magnetic field of 
0.1 T applied parallel to the Cu02 planes where the indium 
counter electrode is in its normal state without field inducing 
spectral peaks. The red curve is for 1 T in increasing fields 
applied perpendicular to the Cu02 planes, where field induced 
nodal energy scale peaks are observed. 
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FIG. 7: (color online) Spectral peak positions, S, as a function 
of magnetic field for various doping levels films at log- log scale. 
All data shown here was taken in decreasing magnetic field. In 
overdoped films the tunneling spectrum exhibits peaks at zero 
magnetic field, where for underdoped ones peaks are missing 
at low magnetic fields. All 5 values coincide at high magnetic 
fields and follow the dashed line having a slope of 1 meV/H^ . 



pie is slightly underdoped at T c = 88.1K. The black solid 
line shows a sharp zero bias peak without any observable 
splitting. An upper limit for the bias of possible spon- 
taneous peaks is extracted from this measurement to be 
2ksT = 0.06meV '. We note that even at 16 T, we do not 
observe any zero bias peak splitting under this configu- 
ration, which ensures the sample's orientation. By con- 
trast, when the field is applied parallel to the sample's 
c-axis, the two spectral peaks at ±2.2 meV are clearly 
seen. 

In Fig. [5] we show a typical field dependence of the 
nodal energy scale. At low fields up to 1 Tesla, a single 



peak is observed at zero bias. Upon increasing the mag- 
netic field a field induced nodal energy scale appears at 
lower energies when compared to the case of optimum 
doping, until it reaches a cross-over field (marked by H* 
in Fig. [6]) at which it recovers the optimally-doped yH 
behavior. In unserdoped samples there appears to be a 
well defined field below which the ZBCP does not split. 
This field increases rapidly with underdoping. 
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2. Overdoped case 

All overdoped samples exhibit spontaneous zero field 
spectral peaks. Dagan and Deutscher— have shown a 
correlation between the spontaneous spectral peak bias 
values, So, and the rate at which this bias increases with 
field. They concluded that the spontaneous peaks are 
due to a modification in the order-parameter symmetry 
near the surface in the vicinity of optimum doping from 
pure d — wave for underdoped samples, to d ± id xy or 
d ± is in overdoped ones. 

In this section we present a study of overdoped sam- 
ples with different oxygen doping levels in high decreas- 
ing magnetic fields starting from fields as high as 32.4 
T. The results are summarized in Fig. [7J For slightly 
overdoped films (see all data points above the dashed 
line), we find zero field spectral peaks, with a minute 
shift at low magnetic fields. The zero and low magnetic 
field data are qualitatively in agrement with the measure- 
ments of Dagan and Deutscher— . However, new behavior 
is observed at high magnetic fields. At these high fields, 
all data points collapse to a single line having a slope 
of 1 meV/Ha (dashed line). This is the same slope as 
found at high fields for optimally doped and underdoped 
samples. 

It has been suggested that finite energy peaks can re- 
sult from trapped vortices and their associated Doppler 
shifting super-current at the surface^. Another expla- 
nation could be a minority imaginary component of the 
superconducting order-paramete r 15 i 31 . Both cases break 
time reversal symmetry, as the spontaneous currents flow 
in a specific direction. Applying additional currents 
should then result in either increasing or decreasing the 
net current, assuming that the time reversal symmetry 
is broken macroscopically. However, we find no differ- 
ences in the tunneling spectra for both polarities. This 
is shown in Fig. [HI When the magnetic field is increased 
for a zero field cooled sample the spontaneous peak value 
should reduce for about one half of the samples; or when 
the induced current is opposite to the spontaneous one. 
However, in over 100 junctions measured in this study 
we never observed that the spontaneous peak's bias de- 
creased with increasing magnetic field. 

Theoretical studies by Asano et aZi 32 ' 33 and Kalenkov 
et alM- show that surface impurities could also result in 
spontaneous spectral peaks. These models predict that 
scattering of impurities cause bound states that dom- 
inate the low bias spectrum. However the Andreev - 
Saint- James bound states will still be present, and a 
three peaked structure at zero-bias should be present. 
We did not observe such behavior in any of the junc- 
tions presented here. Moreover, according to these mod- 
els, the "impurity" spectral peaks bias value should be 
proportional to the amount of impurities at the surface. 
In the case of a clean junction, a zero bias peak should 
be present. However, when a magnetic field is present, 
a splitting via an Aharonov-Bohm like phase shift will 
occur—. This means that at high bias the zero field 
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FIG. 8: (color online) Spectral peak position, 5, as a function 
of the absolute magnetic field in both polarities (sample S5) . 
We found no evidence that the polarity of the magnetic field 
influences the peak position for a given field. 

spectral peaks will be shifted by an external magnetic 
field and for a clean interface, no magnetic field shifted 
peaks should be observed. However, an earlier study by 
Dagan and Deutscher— (see also Fig. [7J shows the oppo- 
site trend. At low magnetic fields, junctions showing a 
spontaneous spectral peak shift to a lesser degree than 
those showing a zero bias conductance peak. This rules 
out impurities as a possible explanation for spontaneous 
peaks. 



V. DISCUSSION 

As mentioned before both theories - Laughlin's 3 
id xy theory and Landau-states by Gor'kov-Schrieffer - 
Anderson^ result in exactly the same field dependence 
for the induced nodal scale. However, there are two 
main differences between these approaches that can be 
checked experimentally. First, Laughlin predicts a weak 
first-order phase transition to the id xy state which is 
not predicted by the Landau-state theorem. This phase 
transition was in fact demonstrated recently by Elhalcl 
et al^. Second, the Gor'kov-Schrieffer - Anderson the- 
ory predicts a series of states while in Laughlin's theory, 
only one energy scale appears. The second peak is not 
observed, even down to 20mK. Additionally, we showed 
that the trajectories between two successive Andreev - 
Saint- James reflections are much longer than the films 
thickness. It is therefore unlikely that such states exist 
in the thin films used in our measurements. 

Laughlin's theory however, has no doping dependence, 
which is a key feature observed in our measurements. Fol- 
lowing Laughlin, Deutscher et. alM suggested a doping 
dependence correction to the free energy in the form 
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F = aS 2 + b8 3 - cSB (3) 

Here b and c are calculated by Laughlin. a is a doping 
dependent term, a — cio(x — x c ), where ao is a negative 
constant and x c is the optimal carrier concentration. 

Using Eq. [3l we calculate a minimum F for (5^0 
at zero field only in the overdoped regime where x > x c , 
while at higher fields the square root behaviour of S is re- 
covered for both underdoped and overdoped regimes. We 
can therefore conclude that our data is better described 
by the modified Laughlin's theory 3 -^ with an additional 
order parameter component. 

We have shown that time-reversal symmetry is not 
broken, macroscopically, even when spontaneous spec- 
tral peaks appear in the tunneling measurements. An 
experiment, similar in concept, was conducted by Tsuei 
et al. where they measured the spontaneous half flux- 
quantum vortex in a tri-crystal experiment and found 
no difference between spontaneous vortices at opposite 
polarities^. The tri-crystal experiment claimed to rule 
out a minority component to the superconductor order- 
parameter. 

Because our measurements, as well as the tri-crystal 
ones, are macroscopic, domain regions with alternating 
spontaneous current directions can reconcile both experi- 
mental results. The origin of such regions could be alter- 
nating minority order-parameters having ±id xy or ±is 
symmetries as both plus and minus states are degener- 
ate. In fact, such configurations could be energetically 
favorable. In such a case, at the domain wall region the 
spontaneous current should be zero. Our technique has 
the advantage that it can probe the zero magnetic field 
state, and the microscopic time-reversal symmetry break- 
ing, regardless of the spontaneous currents direction at 
the microscopic scale. 

Since the order-parameter changes over a length scale 
set by the superconductor coherence length, one should 
be able to find nanometer scale regions (the domain wall 
region) where the minority component order-parameter 
is zero, while, in other regions (inside the domains), it 
should be finite. Only a microscopic study, for example 
with a scanning tunneling microscope, may be able to 
probe such domains. Furthermore, in the domain wall 
region, the order-parameter symmetry should be purely 
e?-wave. Therefore, measurements aiming at detecting 
node- line excitations, such as thermal conductivity, may 
be dominated by the domain wall regions. 



VI. SUMMARY 

In conclusion, tunneling experiments revealed that the 
spectrum of quasi-particle states in nodal regions of a 
(i-wave superconductor is profoundly modified by apply- 
ing a magnetic field perpendicular to the Cu02 planes. 
Doppler-shift of the field induced nodal energy scale has 
been identified, and as expected, is large enough to pre- 
vent observation of a field induced nodal energy scale, 
however it has been minimized by choosing an appropri- 
ate geometry. We showed that the zero field spectral 
peaks cannot be explained by either inelastic scatter- 
ing or trapped vortices at the surface, but rather by a 
domain-like structure of minority order-parameter com- 
ponents with alternating signs. We studied the interplay 
between the spontaneous spectral peaks with the forma- 
tion of the field induced nodal energy scale and film dop- 
ing. Although the low energy states are in agreement 
with theories based either on the explicit description of 
Andreev - Saint- James reflections by the order parame- 
ter away from the nodes^, or on a free-energy expres- 
sion that takes into account a gain in energy due to the 
interaction between the applied field and the magnetic 
moment created by an id xy component 3 . But the ab- 
sence of higher energy level peaks and unreasonably long 
length of the trajectories that would be necessary to ob- 
serve Landau levels are in contradiction with the former, 
while the doping dependance and the first order phase 
transition observed by Elhalel et aL— , are in favor of the 
existence of a field induced id xy order parameter. 
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